Power spectrum of electron number density perturbations 
at cosmological recombination epoch 



Venhlovska B., Novosyadlyj B. 
February 19, 2009 

O : 

^ ' Astronomical Observatory of Ivan Franko National University of Lviv 

X5 : 

<D ■ Abstract 

^: 

, The power spectrum of number density perturbations of free electrons is obtained for 

^—t ■ the epoch of cosmological recombination of hydrogen. It is shown that amplitude of the 

■ electron perturbations power spectrum of scales larger than acoustic horizon exceeds by 
^ I factor of 17 the amplitude of baryon matter density ones (atoms and ions of hydrogen and 

S^' helium). In the range of the first and second acoustic peaks such relation is 18, in the 

O ' range of the third one 16. The dependence of such relations on cosmological parameters 

^ ! is analysed too. 

: Introduction 

> : 

^ ■ The observational data on the cosmic microwave background (CMB) temperature fluctuations 
! obtained in the ground-based, balloon and space experiments became the key test of cosmolog- 
^ I ical models of the Universe and the main source for the determination of their parameters. The 
CN I data obtained in the space mission WMAP (Wilkinson Microwave Anisotropy Probe) lU El [3] 

■ are the most outstanding achievement of the modern cosmology. They have given us the pos- 
O ■ sibility to estimate the cosmological parameters (see [H E] and citing therein) with accuracy 

>■ , that practically equals to one of theoretical predictions because of the non completeness of 
^ ! physical processes governing the CMB temperature fluctuations and polarization as well as the 
^ [ accuracy of analytical approximations and numerical computations. The important part of the 
theory of CMB anisotropy is the kinetics of cosmological recombination of primary plasma and 
calculation of free electrons number density at decoupling of thermal electromagnetic radiation 
and matter. The bases of the theory of cosmological recombination have been founded by 
Zeldovich |7| and Peebles [8] in 1968. In the following papers O [IHl [IH [H, [HI [H] the main 
processes were particularly studied and a few-percent accuracy of calculations of recombination 
kinetics was achieved. The most complete analysis of the kinetics of cosmological recombination 
was carried out by Seager et al. (2000) [I5], where the multi-level (~ 300 levels) structure of 
hydrogen and helium atoms, non-equilibrium kinetics and, practically, all processes determin- 
ing the thermal state of plasma were taken into account. The authors have created also the 
publicly available code RECFAST |16f1 . which provides the computation of electron number 
density with accuracy ~ 1%. It is used in the CMBFAST ^ [18], CAMBCODE p], CM- 
BEASY [20] codes for calculation of power spectrum of matter density perturbations as well as 



-'the improved last release was done in September 2008 
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CMB anisotropy. Since the accuracy of the WMAP observational data is close to the theoretical 
predictions one and is expected to be few times higher in the planned experiment PLANCK, 
the analysis of physical processes not taken into account until now becomes actual. Among 
them are the transitions between the high levels of atoms of hydrogen and helium, their fine 
structure and elaboration of spontaneous, collision and induced transition rates, absorption and 
scattering of thermal radiation along the line of sight, isotope composition, molecular formation 
and dissociation, etc. In recent years these and other effects have been actively studied in order 
to improve an accuracy of computation of cosmological recombination [211 ESI ESI EH [25]. In 
our papers [26l [27] the influence of adiabatic density perturbations on the number densities of 
ionized fractions has been studied. It was shown that at decoupling epoch the amplitude of 
relative perturbations of electron number density is 4-5 times higher than amplitude of relative 
perturbations of total baryon density. The reason for such difference is the considerably higher 
sensitivity of hydrogen and helium photoionization rate to temperature fluctuations than the 
spontaneous recombination rate one. It is clear, since the ionization of hydrogen and helium 
is provided by quanta of background thermal radiation of Wien spectrum range {hu 3> kT). 
The higher amplitude of relative density perturbations of free electron component, meanwhile, 
is not imprinted in the amplitude of angular power spectrum of CMB temperature fluctuations 
in the linear approach^. However, it is not excluded that this effect has imprint in the map 
of the CMB primary polarization. Therefore, in this paper we study more detailed this effect 
for adiabatic perturbations of different scale at region of acoustic peaks to compare the power 
spectra of number density perturbations of different compoments at decoupling. All computa- 
tions were carried out using the code drecfast.f which is described particularly in [26], [27] and 
is publicly available at \hUp:/ /astro. franko.lviv.ua/^novos/\ Research of the paper is restricted 
by the ACDM-model with parameters determined on the basis of data on CMB temperature 
anisotropy, large-scale structure and dynamics of expansion of the Universe. Combination of 
the WMAP data [2] with different datasets on large-scale structure of the Universe, dynamics 
of its expansion, etc. gives somewhat different values of parameters. The ranges of values of 
cosmological parameters determined by [4] are as follows: = 0.7 ^ 0.8, Vt^ = 0.23 -i- 0.31, 
VLb = 0.04 -^ 0.05, h = 0.68 ^ 0.75, = 0.75 0.92, = 0.9 ^ 0.96, where = A/3i^o^ 
Qh = Pb/ Per = Pm/ Per eiXQ cosmological coustaut, baryou matter density and total mat- 
ter density (baryons + dark matter) in units of critical density per = SHq/StiG respectively, 
h = i^o/lOOkm/s/Mpc is dimensionless Hubble constant, Ag is amplitude of initial power spec- 
trum of matter density perturbations, Ug is spectral index of scalar mode of perturbations. 
Computations of cosmological recombination and power spectra of electron number density 
perturbations will be carried out for ACDM-model with two sets of the best-flt parameters: 
= 0.736, rim = 0.278 fife = 0.05, h = 0.68, as = 0.73, = 0.96 [5] and = 0.76, 
rim = 0.24 rib = 0.042, h = 0.73, A, = 0.83, = 0.958 \4\. The angular power spectra of CMB 
temperature fluctuations calculated for them are practically identical [5]. 

1 Dependence of relative number density of free electrons 
on redshit in the ACDM-model 

Important feature of the cosmological recombination of helium and hydrogen in the ACDM- 
model is that it occurs when the total energy density of thermal radiation approximately is 
comparable to baryon matter one: ~ c^pb- Henceforth we will use the following definitions: 
Tiui and null are number density of neutral and ionized hydrogen atoms; nnei, '^Heii and nneiii 

^for second order effect see, for example, Khatri & Wandelt, 2008, arXiv: 0810.4370 [astro-ph] 
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are number densities of neutral, singly and double ionized helium; Ue = nuu + nueii + 2nHeiii 
is number density of free electrons; nn = ^hi + '^hii is total number density of hydrogen 
nuclei; nne = ^Hei + '^Heii + '^Heiii is total number density of helium nuclei. We use the relative 
number densities (ionization fractions): Xhi = ^hi/^h is relative abundance of neutral hydrogen, 
xhii = nmi/nu are relative abundances of ionized hydrogen, XHei = '^Hei/'^He, a^Heii = '^Heii/'^He 
and XHeiii = nielli /riue- relative abundances of neutral, singly and double ionized helium, 
Xe = Ue/nu - relative number density of electrons. The ratio of total number densities of 
helium and hydrogen nuclei we define as fue = ^He/^H, which can be expressed via mass 
fraction of primordial helium Yp, so that fue = ^p/4(l— Ip) (further we assume Yp = 0.24 |28]). 
These quantities obey obvious relationships: Xe = xmi + fueXueU + 2/Hea^Heiii, xm + xmi = 1, 
XHei + a;Heii + s^Heiii = 1- The total mass density of baryons can be expressed via hydrogen 
number density and mass fraction of primordial helium in the following way: pb ~ mpUb = 
mp{n}i + 4?T,He) = "^p^h(1 + 4/He), where is mean number density of baryons (protons and 
neutrons), where rrip is mass of proton. 

At early stages of evolution of the Universe {z > 10'^) all hydrogen and helium atoms were 
ionized completely by thermal photons, so xhii = 1, xui = 0, XHeiii = 1, x^ei = xueu = 
and Xe = I + 2/He [m [H]. At 2; ~ 8000 thermal photons with energies higher than ionization 
potential of Hell from ground level reside in the short-wave tail of Planck function and their 
number density becomes too low to keep all helium in the ionization state of Helll. It begins 
to recombine and Hell ions appear. Recombination of Hell occurs in the conditions of local 
thermodynamic equilibrium (LTE), so, the ionization fraction of helium, XHeiii, is described by 
Saha equation: 

XeXHeiii ^ (27rmeA;r^)3/^ ^_y^^„/;,r„^ 
XHeii h^nu 

where Tm is matter temperature, rrie is mass of electron, h is Planck constant, k is Boltzmann 
constant, XHeii is ionization potential of Hell. Since at this epoch both hydrogen and helium 
are completely ionized (xm = 0, xhii = 1, a^Hei = 0) we have XHeii = 1 — a^Heiii and Xg = 
1 + /He(l + ^^Heiii), SO cquatiou ^ can be easily solved for Xg. Using it one can easily check 
that already at z ~ 5000 all helium atoms become singly ionized. Such state is kept up to 
z ~ 3500 when Hel begins to recombine. At this time the conditions are close to LTE yet. The 
metastable 2s level plays insignificant role in deviation of radiative recombination rate of Hel 
from LTE one until the part of Hel is less than 1% of total helium content and ionized fraction 
XHeii is described yet enough accurately by Saha equation: 

^e^^Heii ^ ^ (27rmeA:T^)=^/^ ^_^^^^/fer^^ 

where XHei is ionization potential of Hel. Now XHeiii = and x-^ei = 1 ~ ^Heii- For accurate 
calculation of XHeii we must have the exact value of Xe = xhii + f-ReXneii- Despite xhii ~ 1, 
the decrease of riuu in 0.1% caused by the hydrogen recombination leads to the variation of 
He comparable to one caused by Hel recombination because of the domination of hydrogen 
(/h = nii/{nu + riue) = 0.921). So, at this step the hydrogen recombination already must be 
taken into account too. The ionized fraction xhii is described yet enough accurately by Saha 
equation: 

^eXmi ^ (27rmefcr^)3/2 ^_^^^^,^^^ 
xm h^nu 

where xhi is ionization potential of HI. The system of these two equations can be reduced to 
the single cubic equation for Xe, which has one real root: 



Xe = 2J-^/3cos(a/3) -B/3, (4) 
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where B = i?Hi+-RHei, Rnei and -Rhi is right hand of equations (I2|) and ([3]), cos a = C /2^—A'^/27, 
A = D — B'^/3, D = i?Hi-RHei — -Rhi — /He-Rnei, C = 2B^/27—BD/3 — E, E = —RujRneii^ — fue)- 
The code RECFAST was complemented by this solution in order to achieve more accurate com- 
putations of number density perturbations of ions. However, it does not change the results of 
calculations of unperturbed x's noticeably [26l [27]. 

Metastable levels 2s Hel and HI cause a delay of recombination of HeH— >HeI and HH— s>HI, 
violation of LTE population of levels and equilibrium of ionization-recombination processes 
("bottleneck" effect). Saha equation does not already describe adequately the recombination 
and equations of detailed balance must be used [16] : 



and 
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"Hel = q 




1-p 




m^c \ 



(7) 



au = F ■10-'''af/{l + cf),m:'c-' t = T^/W. (8) 



OHei and ctH are effective recombination coefficients of helium [29] and hydrogen [30], respec- 
tively. -fCHei = -^Hei2ip/[S^-^(-^)] ^^c factor taking into account the cosmological redshifting 
of Hel 2^p— V-s photons and = A^2p/[8vrif(2;)] is the factor taking into account the cosmo- 
logical redshifting of Lya photons. The effective photoionization coefficients in ([5]) and ([6]) are 
calculated via effective recombination coefficients as follows: 

p = a{27cmekT^/hy/^e-'"'^'-'^/''^'-. (9) 

Before z ~ 800 the matter temperature Tm practically equals radiation temperature Tr 
since until the time-scale of Thomson scattering remains essentially lower than the time-scale 
of expansion of the Universe, tT/feubbie < 10~^. Therefore, the rate of temperature decreasing 
is governed by adiabatic cooling of radiation (7 = 4/3) caused by expansion of the Universe: 

"d7 " (TTI) ■ ^ ' 

And only after recombination at z < 800 adiabatic cooling of ideal gas (7 = 5/3) begins 
to dominate over the heating caused by Compton effect which is the main process of energy 
transfer between electrons and photons. Cooling of plasma via free-free, free-bound and bound- 
bound transitions and collisional ionization as well as heating via photoionization and collisional 
recombination gives insignificant contribution into the rate of temperature change, it does not 
exceed the 0.01% of main processes - adiabatic cooling and heating by Compton effect [T5] . 
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Figure 1: Top panel: the dependence of relative number density of free electrons on redshift in 
ACDM-model with parameters = 0.736, = 0.278, Vlb = 0.05, h = 0.68 f5] (solid line) 
and VIa = 0.76, D.^ = 0.24, fi^ = 0.042, h = 0.73 [4J (dashed line overlaps with solid one). 
Dotted line represents visibility function dr/dze''^ (x270). The relative difference Xe for two 
sets of parameters of ACDM-model is shown in bottom panel. 



So, at this epoch the following equation for rate of temperature decreasing is enough accurate 

dz ~ 3H{z){l + z)meCl + Me + Xe^ (1 + ^) ' ^ ' 

The Tabld2]lists the values of all atomic constants and coefficients used in the equations (IB- lfTTIl . 

In FiglHthe relative number density of electrons Xe at range 400 < z < 10000 is presented. 
The visibility function dr/dze''^ (dotted line) shows the region of the largest cosmological 
recombination rate and decoupling epoch - here r is optical depth caused by Thomson scattering 
by electrons, z = (a~^ — 1) is redshift. 

The calculations shown in FigH] have been done for ACDM model with parameters f^A = 
0.736, rim = 0.278, fif, = 0.05, h = 0.68 [5] and VIa = 0.76, = 0.24, Vlb = 0.042, h = 0.73 [4]. 
The curves Xe{a) for two sets of parameters practically overlap: the difference ~ 4% is in the 
range of maximum of visibility function {zdec ~ 1080) and increases up to ~ 6% for residual 
ionization Xe ~ 10^^ at 2; < 900. 



2 Perturbations of number density of ions and electrons 
2.1 Definitions 

Let the values of ionization fractions of HI, HII, Hel, Hell, Helll and electrons averaged over the 
whole space at fixed cosmological time be Xi, where "i" marks each component. Let us denote 
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the local value of relative number density of each component in the range of cosmological 
density perturbation of baryon matter 6b = Spb/pb -C 1, were pb is mean matter density of 
baryons, as Xj. Its deviation from mean value we denote by 6xi, so that Xi = Xi + 6xi and 
6xi is called the perturbation of relative number density (fraction) of the i-th component. 
We define relative perturbations of ions and free electrons fractions as Aj = 6xi/xi, where 
Xi is any of them. It is obvious that Ag = Sne/ue — ^nn/nn, Ann = Srimi/nmi — Snu/nu, 

Anell = ^'^Hell/'^Hell-'^'T'He/'^He, Anelll = ^?^HeIII/^^•HeIII " ^'^He/'^He- Wc SUppOSC the primordial 

chemical composition of baryon matter to be uniform (/ne is constant) and for homogeneous 
mediunj^ Snu/n-n = Sriuc/nuc = ^b-, then 

Ai = 5, - 6b. (12) 

where 5i = Sni/rii = Ai + Sb is relative number density perturbation of i-th component. It 
must be noted that in expanding Universe the recombination does not end with the completely 
neutral hydrogen or helium but with residual ionization. Therefore, none of values rii reach 
zero and ambiguity of "0/0" -type in Si does not appear. So, Aj's never diverge. Numerical 
results presented in [26l [27| and below prove that. 

Therefore, Aj is difference of number densities relative perturbations of "^"-th component 
and of all baryon matter. Since 6i and Sb are scalar functions of four coordinates in some 
gauge, under the gauge transformations not changing the cosmological background each of 
them is transformed by adding the same expression from the time coordinate transformation 
component (see, for example, [HH [32l [33]). As soon as they appear in ( fT2ll with opposite signs, 
Aj's keep unchanged under such transformations, so they are gauge-invariant variables. 

If hydrogen and helium are entirely ionized and ionization degree does not change with time 
then 6i = 6b and Aj = 0. If the photorecombination and photoionization rates as well as ioniza- 
tion degree of some component change in space and time then 6i and 6b can evolve with different 
rates because the variation of 6b is driven by gravitation and stress of baryon-photon plasma 
and 6i is additionally influenced by kinetics of ionization-recombination processes. Therefore, 
Ai is a measure of deviation of relative number density perturbation of "«"-th component from 
relative density perturbation of total baryon component 6b, caused by different recombination 
and ionization rates within region of cosmological density perturbation. 

At early stage of the Universe evolution the adiabatic relative density perturbations of 
baryon matter 6b and radiation energy 5r = 6e^/eii obey the following relation: eR = 45^/3. 
Since eR = aT^, 6Ta = STr/Tr = 1/35;,. 

2.2 Equations 

The amplitudes of the baryon density 6b and thermal radiation 5^^^ perturbations generated 
in the early Universe increase because of gravitational instability and at the moment of re- 
combination they achieve a value ~ 10~^ — 10~^ at scales 30-300/i~^Mpc (they also depend on 
power spectrum of initial perturbations). The local baryon mass density perturbations lead 
most probably to corresponding perturbations of number density of ions and electrons, 6e oc 6b. 
Since the rates of ionization-recombination processes depend on density and temperature of 
baryon matter and radiation, in the region of the perturbations the departure of the distribu- 
tion of atoms over ionization states from the background one will occur, so A, 7^ is expected. 
We study the cosmological perturbations of small amplitudes. It means that within region 
of cosmological perturbations all equations (fTl)- (fTTI) are applicable and connection between the 

^It is also macroscopically electroneutral: Ue — np + tihcII + 2nHciii everywhere and always. 
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perturbations of ion number density and cosmological perturbations of density and temperature 
can be obtained by variation of those equations. 

Varying the equation for electron ionization fraction Ae = XHeiii/HeAHeiii/a^e we got: 



^ _ a;Heiii(l — a;HeIIl)/He 



Xe + {I- a;Heiii)a;Hciii/Hc 



3 + XHcn\ _ 

2 kTm) 



(13) 



We see that relative perturbation of electron number density is linear combination of initial 
relative perturbations of baryon matter temperature and density. In the region of adiabatic 
perturbations the fluctuations Ag and 6t,^ have the same sign and opposite one to the baryon 
density perturbation 6b. The values of Xe and XHeiii are calculated from ([T]). The asymptotical 
behaviour of Ag follows from (fT3l) : at 2; > 7000 when XHeiii 1 (all helium atoms become 
double ionized) Ahciii — ^ ((5nHoin — ^b) and at redshift z < 5000 when XHeiii (all helium 
atoms become singly ionized) Ag — > 0. So, Ag has peak in this range of redshifts (see Fig. 2 in 

m)- 

At 3500 < z < 5000 both hydrogen and helium are entirely ionized (helium singly): xmi = 
x-Reii = 1, = xuei = Xueiii = 0, SO the amplitudes of all relative perturbations equal to 
zero. With subsequent decreasing of temperature Hel atoms and afterwards HI ones begin to 
recombine. The kinetics of their recombination is described by Saha equations ([2]) and ([3]). 
Variation of these equations gives the expressions for relative perturbations of helium Aneii 
and hydrogen Ann fractions, using them the relative perturbation of free electrons fraction Ag 
can be presented in the form: 

(1 - xmi)xnu^ (1 + (1 - XHeii)a;Heii/a;e) 

Ae = —J- ^ ^ 7. (>T^ + 



+ 



1 - xmi)xmi + (1 - a;Heii)a;Heii/He - Xe 
[1 - XHeii)a;Heii|^ (/hb - (1 " a;Hii)a;Hii/a;e 

(1 — a;Hii)a;Hii + (1 — a;Heii)a;Heii/He — x^ 
[1 - a;Hii)a:Hii + (1 - xucn)xRenfHe 



1 — XmijXiiu + (1 — a;HeIl)a;HeIl/He — Xe 



—6t — Si 



(14) 



At XHeii 1 and xmi — > 1 Ag as expected. Another asymptotical behaviour (xHeii and 
a^Hii 0) has not physical sense as soon as at xneii < 0.99 and xhii < 0.99 it is necessary to use 
the non-equilibrium kinetics equations and energy balance ([5|1- (|TT1) . In this case the differential 
equations for relative perturbations Ann, Aneii and can be obtained by the variation of 
(l5!l-(fTTi). Such equations in the explicit form are presented in [26l[27], their generalized form is 
following: 

= ^ + {t^ {B, -Q + A) +A-l}\ + {A, + Q- A} Sb+ 

(i - - + ^^A&) (i + fef) - B^'-m + e^fe - A)} 5tJ . (15) 

where index of "i" has two values corresponding to Hell or HII and coefficients Ai, Bi, Ci, Di, 
Fi, Bj are determined by number density of ions of helium or hydrogen and the recombination 
rates (see Appendix). The equations of non-equilibrium recombination ([5]) and (l6|) are written 
as follows: 

dxj _ XiXeUnai Dj Aj 

dz ~ H (z) (1 + z) AiQ A, + A ■ ^ ' 

The expression for relative perturbations of free electrons fraction can be found from (fT5l)-(fT6l): 

Ae = ^S^^gf^-i-[{T^(i?.-a + A)+A.-l}A. + {A. + a-A}5.+ 
{f^ (l - Ax5s) - + D^A) (1 + Sf) - ^^ftf + 0* - A)} 5tJ . (17) 
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The evolution of relative perturbations of baryon matter temperature is described by equation: 

Tm'-^ = sHt)^%Li^f:U.. [(Tm - Tk) t^^A, + (4T^ - 5Tk) 5r, + Tr5^J . (18) 

Thus, the system of first-order ordinary linear differential equations for relative perturba- 
tions of matter temperature and ions and electrons fractions consists of equations (fTSll-ffTSi) and 
can be solved using the publicly available code DVERK [34] . 

The equations ffT3l)- (fT8l) were used to analyse the evolution of relative density perturbations 
of ions and temperature perturbations of baryonic matter. All equations contain the solutions of 
unperturbed problem therefore it seems naturally to supplement the code RECFAST [16] with 
block calculating the perturbations of ionization fractions. The complemented code drecfast.f 
[39] is used further in our analysis of evolution of number density perturbations of free electrons. 



3 Evolution of relative density perturbation of free elec- 
trons 

Before recombination the time variations of baryons density and thermal radiation temperature 
perturbations depend on the relation of scale of perturbations to acoustic horizon scale [35] . 
When the scale of perturbation becomes substantially smaller than scale of acoustic horizon 
(Jeans scale) before recombination, then adiabatic perturbations in the baryon-photon plasma 
start to oscillate like the standing acoustic waves. In consequence of recombination the Jeans 
scale drops and the previously oscillating amplitudes of perturbations in baryon component 
start to increase monotonously mainly as a result of gravitational attraction of dark matter 
density perturbations. The amplitudes of perturbations with scales larger than acoustic horizon 
at recombination epoch increased as 5b oc t^^'^ in radiation-dominated epoch and 5b oc t^^^ 
after recombination in dust-like Universe. In the papers [3ll [32], [36l [37l |33l [38| one can find 
the analytical solutions of relevant equations for evolution of relative density perturbations in 
simplified cases of single component media as well as the numerical solutions for real multi- 
component Universe. 

Let us calculate the evolution of number density perturbations of free electrons 5e in the 
region of positive initial matter density perturbation {5b{zinit, k) > 0, Zinu 2> Zdec) of different 
scales in the ACDM-model for the range of redshifts 200 < z < 10000. For this we will integrate 
the system of equations (fT3l) -f|T8ll using the code drecfast.f . Results of calculations of number 
density perturbations of free electrons 5e, baryon density perturbations 5b and radiation density 
perturbations for adiabatic perturbations with wave numbers k = 0.01, ...0.1 Mpc~^ are shown 
in FiglSl The visibility function dr/dze''^ is also shown, its peak denotes the position of last 
scattering surface. In ACDM-model with parameters [5] it is at Zdec = 1080, the wave number 
corresponding to radius of acoustic horizon at this moment equals k^g^ ~ 0.037Mpc~^, to 
particle horizon is kdec — 0.021Mpc^^. 

As one can see in FiglH before beginning of recombination of hydrogen {z > 1500) the 
amplitudes and phases of relative perturbations of free electrons 5e and baryons 5b number 
density coincide for all scales: Ag ~ 0. At the epoch of cosmological recombination the relation 
between them depends on a scale of perturbations. If k < k^g^ then amplitude of electron num- 
ber density relative perturbations is approximately 4 times higher than amplitude of baryon 
matter relative perturbations. For smaller scales, k > k^gg, difference between the values of 
amplitudes of free electrons and baryons relative perturbations is determined by phase of oscil- 
lation of temperature perturbation. After recombination 5e < 5b for all scales. It is so because 
we analyse here the adiabatic positive initial perturbations for which the cold dark matter den- 
sity perturbation increases all the time and after recombination baryon matter falls into the 
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Figure 2: The evolution of number density perturbations of free electrons in region of adiabatic 
matter density perturbations in ACDM-model [5] for k = 0.01, ...O.lMpc"^. (The figure for 
other set of k one can download from |http:/ /astro.franko.lviv.ua/^novos/fig2.pdf| ) 
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Figure 3: a) The relation of amplitude of free electrons number density perturbations and 
amplitude of perturbations of total baryon density, 5e/5b (solid line), the relation of amplitude 
of electron density perturbations and amplitude of radiation perturbations, 5e/^R (dotted line), 
the relation of amplitude of perturbations of total number density of baryons and amplitude 
of radiation perturbations, 5h/5R (dashed line), b) The relation of amplitude of free electrons 
number density perturbations and amplitude of baryon number density perturbations 5e/5b for 
ACDM-model with parameters = 0.736, fi^ = 0.278 fif, = 0.05, h = 0.68 [5] (solid line) 
and = 0.76, Vt^ = 0.24 n^, = 0.042, h = 0.73 [4] (dashed line). The relative difference in 
percents is shown by dotted line. 

potential well caused by cold dark matter perturbations. So, shortly after recombination the 
values of baryon density perturbations acquire the same sign and values that CDM ones and 
practically do not depend on phase of oscillation at the moment of decoupling. 

The main part of the CMB photons was scattered by free electrons in the region of maximum 
of visibility function. Thus, the relation of amplitudes of relative perturbations of free electrons 
and baryons number density, 6e/Sb, at Zdec can define some features of CMB anisotropy. That's 
why a more detailed analysis will be made exactly for this time moment. In FigjSh the relation 
of amplitudes of electron and total baryon number density perturbations (5e/5b) is shown by 
solid line, the relation of amplitudes of electron and radiation number density perturbations 
(Se/Su) by dotted one and the relation of amplitudes of baryons total number density and radia- 
tion perturbations {Sb/Sn) - by dashed one for range of scales 0.001 < k < 0.15. For scale larger 
than particle horizon (k < k^ec), the relations are approximately scale-independent: Se/Sb ~ 4.2, 
^e/^R ~ 3.2 and Sb/S^ ~ 3/4. At lower scales (k > kdec) they change in wide ranges - this is de- 
termined by the different oscillation phases of electron, baryon and photon components pertur- 
bations at Zdec for different scales. At Figl3b the peaks correspond to close to zero values of am- 
plitudes 5b and 5r (for 5b zeros are at k ^ 0.0298, 0.0485, 0.0718, 0.0915, 0.114, 0.134 Mpc"^ 
and for (5r at k 0.0296, 0.0483, 0.0710, 0.0907, 0.112, 0.132 Mpc'^). Zeros of 5e are dis- 
placed to lower scales comparing with them: k^ 0.0305, 0.0501, 0.0737, 0.0941, 0.117, 0.137 
Mpc-^ 

To estimate the dependence of the 5e/5b at recombination epoch on values of cosmolog- 
ical parameters we calculated this ratio for two ACDM-models mentioned in Section 1 (see 
also FigH]). The results of calculations are shown in FigjSb. One can see, that for super- 
horizon perturbations the ratio 5e/5b practically do not depend on parameters of cosmological 
model (the difference is less than 0.5%). For lower scales the difference is significant (> 10%) 
and increases with decreasing of scale. At zero point of 5b the "jumps" are caused by their 
displacement. In ACDM-model with parameter [4] zeros of 5b at the peak of visibility func- 
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tion correspond to scales k ^ 0.0302, 0.0489, 0.0726, 0.0924, 0.116, 0.135 Mpc~\ and zeros 
of 5r to scales k ^ 0.0300, 0.0488, 0.0718, 0.0917, 0.114, 0.134 Mpc'^ The zeros 6e are at 
k^ 0.0310, 0.0506, 0.0745, 0.0951, 0.118, 0.139 Mpc'^ 



4 Power spectrum of number density perturbations of free 
electrons 

The calculation of the power spectrum of density perturbations of any components at any time 
moment t requires the calculation of the transfer function which is defined as follows: 

Ti{k,t) = 5i{k,t)/Si{kmin,t), 

where kmin k^ec and 6i(k,tinit) = 5i(kmm, ^mi*)- It means that transfer function is relation of 
amplitudes of perturbations of two scales k^m and k at any time, initial amplitudes of which 
were equal at initial time tinu- If we have the transfer function than the power spectrum of 
"i"-th component at Zdec can be calculated as follows 

where Ag is normalization constant of power spectrum of scalar perturbations, Ug is spectral 
index. Since we analyse here the relation of amplitudes of electron number density and baryon 
density perturbations, the normalization constant can be arbitrary (free normalization). The 
results of calculations of the cosmological perturbations power spectrum of different components 
are presented in FigH] by the dimensionless magnitude Pj(k)k^ at the moment of cosmologi- 
cal recombination. It summarizes the conclusions deduced from calculations of evolution of 
perturbations for different scales (Figl2]): at moment of cosmological recombination the ab- 
solute value of amplitude of free electrons density relative perturbations is few times higher 
than amplitude of total baryon matter density relative perturbations. It gives also possibil- 
ity for more detailed analysis of the dependence of the both spectra amplitudes relation on 
scale of perturbations. For super-horizon perturbations (k <ti kdec) the tilts of spectra are the 
same for all components Pi(k) oc k"" and relations of their magnitudes are practically con- 
stant: Pfc : Pc : -Pr : -Pe ~ 1 : 1 : 1-8 : 17. At lower scales the power spectrum magnitude 
of photon-baryon plasma oscillates while one of collisionless component (cold dark matter) in- 
creases monotonously. At these scales the relations of power spectrum magnitudes of different 
components depend on scale of perturbations. The positions of maxima in perturbation power 
spectra of photon-baryon plasma components approximately coincide: 

• for thermal radiation - k ^ 0.0175, 0.04, 0.06, 0.08, 0.1, 0.12, 0.1425 Mpc'^; 

• for baryon matter - k ^ 0.0175, 0.04, 0.06, 0.0825, 0.1025, 0.125, 0.145 Mpc'^ 

• for free electrons - k ^ 0.0175, 0.04, 0.0625, 0.0825, 0.105, 0.1275, 0.1475 Mpc"^ 
For positive cosmological density perturbations {Sc{k, Unit) = 5b(k, t 

init) ^ 0) the first maximum 

is due to the perturbations which were in the phase of first maximal compression at Zdec] the sec- 
ond one - to the perturbations which at Zdec were in the phase of first maximal decompression; 
the third one - to perturbations which at Zdec were in the phase of second maximal compression, 
etc. For adiabatic perturbations with (5c(k, ti„jt) = ^^(k, tj„j() < 0) the maxima have opposite 
character - maximal decompression, maximal compression, etc. The positions of "rifts" corre- 
spond to zeros of perturbations magnitude and coincide with positions of corresponding "jumps" 
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Figure 4: The power spectra of number density perturbations of free electrons (solid line), 
baryon density perturbations (dash-dotted), perturbations of energy density of thermal radia- 
tion (dotted line) at the moment of cosmological recombination Zdec- 

in FiglH Since experimental determination of positions and amplitudes of acoustic peaks of 
the CMB temperature fluctuations power spectra and their theoretical interpretation are an 
extremely actual tasks for modern cosmology, we present the ratios of amplitudes of power 
spectra peaks of different components for two ACDM-models in the Tabl^ The comparison 
of the positions and amplitudes of acoustic peaks of density perturbations power spectra of 
all components and zeros of (5r, 6b and 6e for two models shows that differences do not exceed 
~ 1 — 2%. It explains good coincidence of the predicted power spectra of the CMB temperature 
fluctuations with observational data: xLin = 37.2 for ACDM-model with parameters from [5] 
and 37.8 for ACDM-model with parameters from [4] for 33 degrees of freedom of the system. 

For the flrst two acoustic peaks the amplitudes of electron number density perturbation 
power spectrum are ~ 18 times higher than amplitudes of perturbations of baryon matter 
density, for the third one it is ~ 16 times higher. For the next peaks such ratios decrease. 

Here we do not discuss the relations of amplitudes of baryon, cold dark matter and thermal 
radiation densities perturbations power spectra because they have been discussed in numerous 
works (see, for example [37] and references therein) and are presented here for completeness of 
the picture. 

Conclusions 

At cosmological recombination epoch the amplitude of relative perturbations of electron num- 
ber density, 6e = drie/rie does not equal the amplitude of relative perturbations of baryon 
matter density 5b = 5nb/nb because the recombination and the photoionization rates have 
different dependences on density and temperature of baryon-photon plasma. The difference 
between them becomes prominent when free electrons promptly disappear because of fast re- 
combination of hydrogen at z < 1600 (FigH] and FiglS]). At decoupling of thermal radiation 
from baryon matter the visibility function dr/dze''^ (r is optical depth due to the Thom- 
son scattering by electrons) has maximum at z^tec ~ 1080 in ACDM-model with parameters 
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Table 1: Ratios of peak amplitudes of baryon matter, cold dark matter, thermal radiation and 
free electrons density perturbations power spectra at Zdec for two ACDM-models: [5] - top line, 
[4] - bottom line. 



Np/p 


k [Mpc-i] 


Pb '■ Pc '■ Pr '■ Pe 


1 


0.0175 
0.0175 


1 : 3.7 : 1.70 : 18.2 
1 : 3.5 : 1.72 : 18.0 


2 


0.0400 
0.0400 


1 
1 


67.4 : 1.65 
69.8 : 1.67 


18.9 
18.3 


3 


0.0625 
0.0625 


1 
1 


79.4 : 1.30 
76.2 : 1.34 


16.2 
15.4 


4 


0.0825 
0.0825 


1 : 203.8 : 1.29 : 13.6 
1 : 262.8 : 1.21 : 16.1 


5 


0.1050 
0.1050 


1 : 231.9 : 0.83 : 11.1 
1 : 221.8 : 0.91 : 10.0 


6 


0.1275 
0.1275 


1 : 778.7 : 0.63 : 12.2 
1 : 685.1 : 0.79 : 10.2 



[5]. The power spectra of relative perturbations of electron number density Pe{k,Zdec) = 
{Se(k, Zdec)^eO^, Zdec)) and baryous Pb{k, Zdec) = {Sb{k, Zdec)Sh{K Zdec)) (FigHD computed for 
this moment specify the relation of their amplitudes at different scales. It is practically 
flat for perturbations with super-horizon scales (k <^ kdec)' Pe(X, Zdec) / PbiM, Zdec) ~ 17. On 
sub-horizon scales (k > k^ec) the power spectra of electron number density, baryon den- 
sity and thermal energy density perturbations oscillate. In the ACDM-model the maxima 
of power spectrum magnitudes of electron number density perturbations at moment of de- 
coupling are at scales k^ ^ 0.0175, 0.04, 0.0625, 0.0825, 0.105, 0.1275, 0.1475. For them 
Pe{km, Zdec) / PbiMm, Zdec) ~ 18.2, 18.9, 16.2, 13.6, 11.1, 12.2, correspondingly. Results and 
conclusions weakly depend on change of parameters of ACDM-model. 
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Appendix 

Formulas for calculations of coefficients Ai, Bi, Ci, Di, Fi. 
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A, 



H 12^ 

XYiAiXeUnOiiiei - /^Hei (1 " ^^Heii) e xniiXennoin - /?H (1 - a;Hii) e 



feT„ 



5. 



(3YLe\ (1 - a^Heii) e /3h (1 - Xmi) e" 



XnellXeUnaYlel " /^Hel (1 - a^Hell) 6 '=^™ ' XHIlXgriHaH - /^H (1 " a^ffll) 6 



a. 



-f^HelAHe^H (1 - a;HcIl) 6 iTnAH^H (1 - aJffll) 



1 + i^HelAHeriH (1 " XHell) ' ^ + ^HAnnn (1 - XHIl) 

i^Hel (Ane + /^Hel) (1 - X^iell) 6 fcT™ /^H (Ar + /3h) '^/Z (1 " a;HIl) 



1 + iTHel (Ahc + /^Hel) (1 - XHell) ' 1 + i^H (Ah + /^h) Un (1 - Xffll) 



where f-ReU^s is frequency of Hel 2^s — l^s transition, i/hi2s is frequency of HI 2s — Is transition. 
For helium in contrary to hydrogen it is needed to take into account the slitting of Hel 2^p 
and 2^s that is why there is an additional factor with ^'Hei2ip2is = ^^He^ip ~ ^ndi^s = ^ps in ^ 
comparing to ([6]). The values of function 9j are: Gneii = 1 and 6ofHii = 0. 
Variations of photoionization coefficients were calculated as follows: 

6(3i 6ai 3 hv2si . 

= — Ot H Or , 

Pi Oii 2 kTm 

where hv2si is ionization energy from 2s state. Variations of values of recombination coefficients 
and matter temperature perturbations are connected by 



where 




S —- I -L H / r 



TJT2 1 + Jtjt, J ' V l + c-t\ 
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